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Abstract 

We summarise the perturbative QCD analysis of the structure function data for gi from lon- 
gitudinally polarized deep inelastic scattering from proton, deuteron and neutron targets, with 
particular emphasis on testing sum rules, determining helicity fractions, and extracting the 
strong coupling from both scaling violations and the Bjorken sum rule. 

Talk at Les Rencontres de Physique de la Valle d'Aoste, 
La Thuile, March 1997, and at DIS97, Chicago, April 1997 
to be published in the proceedings 

June 1997 



^ Royal Society University Research Fellow 
^IBERDROLA visiting professor 



The data on the nucleoli structure function gi, determined from measurements of the po- 
larization asymmetry in longitudinally polarized deep inelastic scattering, is now sufficient to 
justify a sophisticated NLO analysis using the full machinery of perturbative QCD It is 
possible to test sum rules, extract polarized parton distributions and determine parton helicity 
fractions, and even, as we shall see, make a reasonably accurate determination of as- Impor- 
tant issues in these analyses are the dependence on the functional forms used to parameterize 
the initial parton distributions, and, when estimating first moments, the implications for the 
extrapolation to small x. Here we will discuss the results obtained in a recent analysis [§. 

In perturbative QCD structure functions are decomposed into convolutions of perturbatively 
calculable coefficient functions and intrinsically nonperturbative parton distribution, which 
then vary with according to perturbative evolution equations |^. For gi the appropriate 
partons are [|, || the nonsinglet and singlet polarized sea quark distributions, Aq^^{x, Q^) and 
Aq^{x,Q^) (where Aq = q^ — q^+q^ — q^), and the polarized gluon distribution Ag{x,Q'^) = g^ — 
gK Under evolution Aq^ and Ag mix, while Aq'^^ evolves independently. Valence distributions 
(of the form q^ — q^ — (g^ — q^)) also evolve independently, but cannot be measured in inclusive 
scattering of virtual photons, and are thus not accessible from data on gi alone (unless of course 
extra assumptions are made about the relative size of valence and sea, as in ref.p, 0). The 
coefficient functions C^^, and Cg and splitting functions P^^, Pgg, Pqg, Pgq and Pgg are 
now all available at NLO P]. 

First moments of polarized parton distribution functions are of special interest, since they 
have a direct interpretation as the helicity fractions of the various constituents. First moments 
of nonsinglet distributions are chosen to be independent of (PCAC), and may then be related 
to nonsinglet axial charges {gA = 0,3 and as for the SU{2) triplet and SU{?>) octet) measured 
in weak baryon decays. The first moment of the singlet distribution AS = /q^ dxAq^ may also 
be chosen to be independent of Q^: this is a direct consequence of the Adler-Bardeen theorem. 
The singlet axial charge is then 

a,{Q^) = A'L-nf^Ag{Q^). (1) 

It evolves multiplicatively, beginning at two loops, because of the axial anomaly: this implies 
that AgiQ"^) grows as logQ^, and that ao(oo) and AS will in general be very different 0. 
Changes in factorization scheme away from these 'AB schemes' result in only small changes 
in Ag (of 0(as)), but large changes in AS (which however is then no longer scale independent). 

Sum rules express the first moments of nonsinglet and singlet combinations of structure 
functions gi in terms of nonsinglet and singlet axial charges: 

Si dxg^ix, Q') = ^[(1 - c^^)(3ea3 + as) + (1 - c^)4ao(oo)], (2) 

where e = -|-1,0, — 1 for N=p,d,n respectively. The isospin triplet (Bjorken) sum rule, for 
example, relates the first moment of gf — g^ to gA- The perturbative corrections c^^ and C^, 
which begin at NLO (i.e. at O (««)), are now both known at NNNLO (i.e. to 0{al)) 

In order to compute first moments of parton distributions and structure functions it is 
necessary to extrapolate the data in the measured region down to x = 0. The behaviour of 
at small x is now well understood due to data from HERA [^: while at low xq^{x, Q^) is 
fairly fiat, in accordance with Regge theory since the pomeron intercept is close to unity, as 
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increases it rises more and more steeply due to a rise in the gluon distribution xg{x, Q^) driven 
by gluon bremsstrahlung |T^. By contrast Regge theory predicts that xq^^{x,Q^) decreases 
rapidly at small x (so q'^^{x,Q'^) -C q^{x,Q^)), though with a rise in q^^{x,Q'^) which is 
then stable under perturbative evolution. Now according to Regge theory |Ag^'^(x, Q^)] and 
\Aq^{x,Q^)\ should be flat (or falling) at small x and Q^, and perturbative evolution then 
leads to an instability in both distributions |13]. However Aq^ and A^f mix nontrivially, and 
the eigenvectors of their evolution are such that they eventually have opposite sign. As a 
consequence although Ag^'^(x, Q^) and A(yf(a:, Q^) both rise at small x and increasing Q"^, 
Aq^{x,Q'^) falls steeply and g^i^x^Q"^) is eventually driven negative 0]. Estimates of small x 
contributions using Regge extrapolation alone will then be unreliable, underestimating their 
size and sometimes even giving them the wrong sign. 

Here our general philosophy will be to estimate small x contributions by performing fits to 
the data which are consistent with Regge behaviour at small x and low Q^, but have the effects 
of perturbative evolution superimposed at higher Q^. Our general procedure is the standard 
one: we take some initial parameterization of Aq^^ {x,Ql), Aq^{x,Ql) and Ag{x,Ql), evolve 
these three distributions to the of the data fl^, |T5|, [1^, |1^ using NLO perturbative evolution 
(in the AB scheme), compute gi{x,Q'^) using the NLO coefficient functions, and then adjust 
the parameters in order to optimize the fit. The quality of such a fit may be judged from the 
figure. In principle it should be possible to fit both triplet and octet nonsinglet distributions 
(the charm and bottom distributions are generated radiatively from threshold), but in practice 
this is not possible with present data since the triplet dominates (og -C 803). We thus assume 
that both distributions have the same shape, and find that the results of the fits are very 
insensitive to the value of ag. 

In such fits it is important to assess the dependence of the results on the form of the initial 
parameterization, otherwise errors may be underestimated. This is particularly true here, where 
we want to compute first moments, since the shape of the small x tails may depend very much 
on the parameterization adopted. To this end we tried four different parameterizations: 

A: the standard form with powers of x and 1 — x, and = IGeV^; 

B: powers of In 1/x and Qg = IGeV^; 

C: as A, but with Qg = O.SGeV^ and falling inputs; 

D: as C, but with valencelike Regge behaviour fixed by hand. 
All four parameterizations give fits to the data of comparable quality. Furthermore they all lead 
to structure functions which eventually become negative at small x (see figure), as expected. 
However the size of the small x contribution varies quite widely between the four distributions, 
and we take this spread as an estimate of the error on the small x extrapolation. It follows 
that our first moments for g^^, gf and gi will all be smaller, but with a larger error, than those 
obtained using Regge extrapolations (see table 1). By contrast our first moment of gi — g^ is 
larger, since the nonsinglet distribution rises at small x. 



Fixing as{mz) = 0.118±0.005 |]T8[, and, very conservatively, ag = 0.6±0.2, we may use the 
fits to determine gA-, AE, Ag and ag. The results are shown in table 2. The theoretical error 
includes the uncertainty in the small x extrapolation (underestimated in ref.[^), and the effect 
of higher order corrections estimated by variation of the renormalization and factorization scales 
(not included in ref. 0). Although the shapes of the distributions are not fixed very precisely 
by the data, their first moments are surprisingly well determined. Fitting gA is a test of the 



2 



Bjorken sum rule: the value obtained from /?-decay is 1.26, so we verify the Bjorken sum rule at 
the 8% level. Similarly fitting AS may be regarded as a test of the constituent quark inspired 
ansatz AS ^ og: again this is confirmed to within one standard deviation (though it must be 
remembered that such a test is only meaningful in AB schemes; in other schemes AS is scale 
dependent and may take any value one chooses). As was predicted in ref. 0, a large positive 
value of A5f(lGeV^) (two standard deviations from zero) is then responsible for the small value 
of the singlet axial charge (compatible results for Ag were also found in P, ^ |6t p!^). 
Our extrapolation to small x is particularly important here: Regge extrapolations give a rather 
larger value for the axial singlet charge, with a smaller (but incorrect) estimate of the error 



1^, There is still no totally convincing explanation of why ao should be consistent with 
zero. 

The recently published neutron data now make it possible to include as in the fit, giving 
the result 

a,(mz) = 0.120;Hg^(exp.);ng^(th.) (3) 

This determination is possible because of the scaling violations in the nonsinglet channel, which 
are not sensitive to uncertainties in the gluon distribution: without the neutron data the 
error would have been twice as large. It is thus not affected by uncertainties in the small 
X extrapolation: the dominant error instead comes from higher order corrections, which are 
particularly large as much of the data is at fairly low Q^. It is tempting to attempt to reduce 
this uncertainty by determining as from the Bjorken sum rule [^, exploiting the NNNLO 
calculations of ref.[|10|. Unfortunately the relatively large error in the Bjorken sum resulting 
from the small x extrapolation makes this method impractical at present: we find as{mz) = 
n 118+0-010 

U.i±0_QQ24. 

The test of the Bjorken sum rule, the determination of helicity fractions, and the result 
for as would all be improved by data with better statistics at higher Q^, by an independent 
measurement of Ag through open charm production [^, or by data at smaller x ^2[. A better 
theoretical understanding of polarized structure functions at small x would also be useful. 
However since the higher order corrections now involve double logarithms of x|^3|], rather than 
the single (and thus factorizable [^) small x logarithms in the unpolarized case, this may 
prove to be a considerable challenge. 
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Regge extrapolation 


pQCD extrapolation 
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0.136 ±0.016 


0.116 ±0.022 
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0.127 ±0.010 


0.107 ±0.017 


SMC 
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10 


0.041 ± 0.007 


0.021 ±0.016 


E143 


d 


3 


0.042 ±0.010 


0.021 ±0.014 


E142 


n 


2 


-0.033 ±0.011 


-0.068 ±0.015 


E154 


n 


5 


-0.041 ±0.007 


-0.070 ±0.015 



Table 1: Results for complete first moments of . 





(jA AE A(/(l) «o('3o) 


Central Value 
Exp error 
Th error 


1.19 0.45 1.6 0.10 
0.05 0.04 0.4 0.05 
0.07 0.08 0.8 tin 



Table 2: Axial charges and helicity fractions. 



5 




Figure 1: Fits to the data: in the left cohimn we plot . and compare fit B to the data, while 
in the right column we plot at = lOGeV^ for fits A-D, showing the differences between 
the small x tails. Figures a,b and c are for proton, deuteron and neutron respectively. 
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